Electromigration in thin tunnel junctions with ferromagnetic/nonmagnetic: 
nanoconstrictions, local heating, and direct and wind forces 



J. Ventur£| and J. B. SousU 



in 
o 
o 

(N 

Cu- 
te 

in 

(N 



c/i 



o 

> 
(N 
l> 

^- 

O 
in 
o 

-I— > 
c3 



i 

c 

o 
o 



X 



IFIMUP and Faculty of Sciences U. Porto, Rua do Campo Alegre,678, 4169-007, Porto, Portugal 

Y. Liu, Z. Zhang and P. P. FreitasOl 

INESC-MN and 1ST, Rua Alves Redol, 9-1, 1000-029 Lisbon, Portugal 

Current Induced Resistance Switching (CIS) was recently observed in thin tunnel junctions with 
ferromagnetic (FM) electrodes i.e FM/I/FM. This effect was attributed to electromigration of metal- 
lic atoms in nanoconstrictions in the insulating barrier (I). Here we study how the CIS effect is 
influenced by a thin non-magnetic (NM) Ta layer, deposited just below the AlOz insulating bar- 
rier in tunnel junctions of the type FM/NM/I/FM (FM=CoFe). Enhanced resistance switching 
occurs with increasing maximum applied current (/max), until a plateau of constant CIS is reached 
for /max ~ 65 mA (CIS~60%) and above. However, such high electrical currents also lead to a 
large (~9%) irreversible resistance decrease, indicating barrier degradation. Anomalous voltage- 
current characteristics with negative derivative were also observed near ± I max and this effect is 
here attributed to heating in the tunnel junction. One observes that the current direction for which 
resistance switches in FM/NM/I/FM (clockwise) is opposite to that of FM/I/FM tunnel junctions 
(anti-clockwise). This effect will be discussed in terms of a competition between the electromigra- 
tion contributions due to the so called direct and wind forces. It will be shown that the direct force 
is likely to dominate electromigration in the Ta (NM) layers, while the wind contribution likely 
dominates in the CoFe (FM) layers. 

PACS numbers: 66.30.Pa, 66.30.Qa, 73.40. Gk, 73.40. Rw, 85.75.Dd 

Keywords: Electromigration, Tunnel Junction, Current Induced Switching, Spin Torque 



INTRODUCTION 



Tunnel junctions (TJ) consisting of two ferromagnetic 
(FM) layers separated by an insulator (1)1 show enor- 
mous potential for a multiplicity of applications such as 
read head^ strain^ current, position and speech sen- 
sors or even to detect magnetically tagged biological 
specimens^ However, probably the most sought after ap- 
plication is high performance, low cost, non- volatile mag- 
netoresistive random access memories (MRAMs) In a 
tunnel junction, the magnetization of one of the FM lay- 
ers (pinned layer) is fixed by an underlying antiferromag- 
netic (AFM) layer. The magnetization of the other FM 
layer (free layer) reverses almost freely when a small mag- 
netic field is applied. Due to spin dependent tunneling! 
one obtains two distinct resistance (R) states correspond- 
ing to pinned and free layer magnetizations parallel (low 
R) or antiparallel (high R). However, several drawbacks 
are still of concern in actual MRAM submicron devices, 
like cross-talk in the array configuration or the large 
power consumption to generate the magnetic field to 
switch R. One then aims to replace the magnetic field- 
driven magnetization reversal by a Current Induced Mag- 
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netization Switching (CIMS) mechanism Such goal 
was recently achieved in magnetic tunnel junction al 11 
for current densities j ~ 10 7 A/cm 2 . On the other hand, 
Liu et alm& observed reversible resistance changes in- 
duced by lower current densities (j ~ 10 6 A/cm 2 ) in thin 
FM/I/FM TJs. These changes, although initially at- 
tributed to the CIMS mechanism, were later founcU^ 
not dependent on the relative orientation of the mag- 
netizations of the free and pinned layers. This effect was 
then called Current Induced Switching (CIS) and is now 
attributecU^*A*i to electromigration (EM) in nanoconstric- 
tions in the insulating barrier. The combination of the 
tunnel magnetoresistive and CIS effects allows the use 
of a magnetic tunnel junction as a three resistance state 
device^ Both CIS and CIMS effects seem to coexist in 
thin magnetic tunnel junctions for j > 10 6 A/cm 2 . The 
reasons for the observed dominance of one effect over the 
other are still unclear but likely related to structural dif- 
ferences in the tunnel junctions. One notes however, that 
electromigration can in fact limit the implementation of 
the spin torque mechanism in actual devices and be a 
major reliability issue in read head sensorsiAl 

When a metal is subjected to an electrical field E, the 
usual random diffusive motion of atoms is biased by the 
resulting driving force F, and a net atomic flux can be ob- 
served. This phenomena is known as electromigrationi& 
and F can be written as: 



Z*eE, 



(1) 
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where Z* is the effective valence and e is the elementary 
charge. The force acting on the migrating ion is usu- 
ally separated into two components, both linear in the 
external applied electrical field: 



(Z d + Z w )eE. 



(2) 



The direct force Fd arises from the electrostatic interac- 
tions between the electrical field and the so called direct 
valence of the ion Zd (> 0). The theoretical calculation 
of the direct force is a challenging process but Zd ~ Z 
(Z =ion valence) is usually assumed. The wind force F w 
results from momentum exchange between the current 
carrying electrons and the migrating ions and so it has 
the direction of the electron current (opposite to the elec- 
trical field) . The wind valence Z w is simply a convenient 
term to describe the wind force, arising from the fact 
that F w is proportional to the current density and, in 
an ohmic material, to E. The competition between wind 
and direct forces is often dominated by the first, which 
usually controls the sign and magnitude of the effective 
valence Z* and the EM process. 

Here we study how a Ta non-magnetic (NM) amor- 
phous thin layer deposited just below the insulating bar- 
rier influences the Current Induced Switching. In a CIS 
cycle, the resistance commutes between two states due 
to electromigration of ions from the electrodes into the 
barrier (decreasing R) and from the barrier back into the 
electrodes (increasing R)pi^ We can then define the CIS 
coefficient as the relative difference between these two 
R-states. Interestingly, the current direction for which 
R-switching occurs in FM/NM/I/FM tunnel junctions is 
opposite to that of FM/I/FM tunnel junctions^ Using 
the intuitive ballistic model of EM, we will show that 
the direct force is likely to dominate electromigration in 
Ta (NM) layers, while the wind force dominates in CoFe 
(FM) layers. The switching direction difference will be 
here associated with the dominance of different EM forces 
(direct or wind) in the two types of tunnel junctions re- 
ferred. 

The CIS coefficient was strongly enhanced by increas- 
ing the maximum applied current (I max ), reaching almost 
60% for /max = 80 mA. However, severe R-degradation 
occurs when I max > 65 mA. Voltage-current character- 
istics show strong anomalous non-linearities, here asso- 
ciated with heating effects. Comparing our experimen- 
tal results with voltage-current characteristics as pre- 
dicted by Simmons' model^ we estimate that the tem- 
perature inside the tunnel junction reaches ~600 K for 
Imax = 80 mA. Numerical results from a model of heat 
generation in tunnel junctions suggest that such high 
temperatures can only occur if local current densities 
much larger than j = I/A (I the electrical current and A 
the total tunnel junction area) exist within the barrier. 
One concludes that these hot-spots concentrate most of 
the current flowing through the tunnel junction stack and 
are likely the reason for the occurrence of EM in the stud- 
ied tunnel junctions. 





FIG. 1: Energy barrier for atomic diffusion, a) without and b) 
with an applied electrical field. Notice how the direction for 
diffusion becomes biased by the driving force F: the energy 
barrier for migration to the right (El) is smaller than that for 
migration to the left {E b ), 



II. ELECTROMIGRATION 

For atomic diffusion to occur, an atom needs to sur- 
mount the energy barrier E b separating neighboring equi- 
librium lattice sites (Fig. ^,). When an electric current 
flows through a metal this usual, thermally-activated, 
random motion of atoms is biased by the electrical field 
(Fig. njj), resulting in a net atomic flow. This phenomena 
is know as electromigrationi& and is currently the major 
cause of failure of interconnects in integrated circuits^ 
Studies of EM in interconnects are performed under 
severe conditions, such as high electrical current den- 
sities (~ 10 7 A/cm 2 ) and temperature (~ 500-700 K) 
and show that EM can occur through different diffusion 
paths, such as grain boundary and interfaces, as in Al21 
and Cu 2 ^ interconnects, respectively. The relative im- 
portance of the different diffusion paths varies with the 
material properties, such as grain size and orientation, 
interface bonding and structure. 

Electromigration is also a concern in magnetic nanos- 
tructures, namely spin valves and tunnel junctions^ 
During device operation, local structural inhomogeneities 
can lead to large current density, and thus to electro- 
migration. This is of particular importance in tunnel 
junctions where the resistance depends exponentially on 
the barrier thickness and where localized nanoconstric- 
tions can concentrate most of the current. Such high 
current densities can also produce intense heating lead- 
ing to enhanced electromigration ■!£ Discrete electromi- 
gration events were observed in metallic nanobridges 
(for j ~ 10 8 A/cm 2 )^ Reversible EM was recently ob- 
served in Ni nanoconstrictions (j ~ 10 13 A/cm 2 )Sl and 
thin tunnel junctions (j ~ 10 6 A/cm 2 )m2^ Electromi- 
gration in these nanostructures can lead to both an in- 
crease and a decrease of the electrical resistance, depend- 
ing on the sense of the applied electrical current, and thus 
on the sense of EM-driven atomic motion. 

The ballistic model of electromigration presents the 
most intuitive picture of the underlying physics of EM. 
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The wind force is calculated assuming that all the mo- 
mentum lost by the scattered electrons is transferred to 
the migrating ioniSi In the free electron approximation 
the wind valence becomes^ 



-nlatr, 



(3) 



where n is the electron density, I is the electron mean 
free path and atr is the electron transport cross section 
for scattering by the ion. Using, e.g. known values for 
Fe (n ~ 1CT 1 A" 3 , I ~ 50 A, a tr ~ 3 A 2 )^. one finds 
Z w ~ — 15 (\Z W \ 3> Z w 2). Such estimative confirms 
that the wind force usually dominates electromigration. 
More elaborated EM models such as the pseudopoten- 
tial method give lower Z w values, by as much as 70% 4& 
However, because of its simplicity, we will use the ballis- 
tic model to qualitatively explain our results. 

Sorbello21 first studied electromigration forces in meso- 
scopic systems. In particular he considered electromigra- 
tion near a point contact, modeled as a circular aperture 
of radius a between two metallic layers of electrical resis- 
tivity p. He found that the direct force is then greatly 
enhanced near such constriction. An estimate on the rel- 
ative magnitude of the wind and direct forces givesiS 3 ^ 



aa tr 

~z7' 



(4) 



which evidences the important role played by the con- 
striction geometry: the smaller the constriction radius, 
the larger will be the direct force compared to the wind 
force. 



The electrical resistance, magnetoresistance and cur- 
rent induced switching were measured with a four-point 
d.c. method, with a current stable to 1:10 6 and using an 
automatic control and data acquisition system. 

CIS cycles were performed using the pulsed current 
methocU^ allowing us to measure the remnant resistance 
of the tunnel junction after each current pulse. Current 
pulses (ip) of 1 s duration and 5 s repetition period are 
applied to the junction, starting with increasing nega- 
tive pulses from I p = (where we define the resistance 
as -Rmitiai), in A / p = 5 mA steps up to a maximum 
+ J maXl dependent on cycle in the 10-80 mA range. One 
then decreases the current pulses (always with the same 
Alp), following the reverse trend through zero current 
pulse (i?haif) down to negative — I ma x, and then again to 
zero (.Rfinai), closing the CIS hysteretic cycle, R = R{I p ). 
Positive current is here defined as flowing from the bot- 
tom to the top lead. 

The junction remnant resistance is measured in the 
5 s-waiting periods between consecutive current pulses, 
using a low current of 0.1 mA, providing a R{I P ) curve 
for each cycle. This low current method allows us to 
systematically discard non-linear I(V) contributions to 
the resistance. However, the voltage across the TJ is also 
measured while applying the current pulse I p , enabling 
us to obtain the (non-linear) V(/ p ) characteristic for each 
CIS cycle. 

Using the definitions above, one can define the CIS 
coefficient as: 



CIS = 



R 



initial 



R 



hall' 



(-^initial + i?half)/2 



(5) 



III. EXPERIMENTAL DETAILS 

In this work we used a series of ion beam deposited tun- 
nel junctions, with a non-magnetic Ta layer inserted just 
below the insulating A\O x barrier. The complete struc- 
ture of the tunnel junctions studied was glass/bottom 
lead/Ta (90 A)/NiFe (50 A)/Mnlr (90 AWCoFe (40 
A)/Ta (20 A)/A10 x (3 A+ 4 A)/CoFe (30 A)/NiFe (40 
A)/Ta (30 A)/TiW(N) (150 A) /top lead. The chosen 
structure is similar to that of magnetic tunnel junctions 
grown for actual applications except for the additional Ta 
layer, thus making a comparison between the FM/I/FM 
and FM/NM/I/FM systems easier. Previous Transmis- 
sion Electron Microscopy images obtained in similar sam- 
ples show no significant microstructural changes induced 
by a Ta layer deposited below the barrier^ The A10 x 
barrier was formed by two-step deposition and natu- 
ral oxidation processes (50 mTorr, 3 min, 100 mTorr, 
20 min)p±£ NiFe, CoFe, Mnlr and TiW(N) stand for 
Ni 80 Fe 2 o, Co 80 Fe 2 o and Mn 78 Ir 22 , Tii W 90 (N). The bot- 
tom and top leads are made of Al 98.5% Si 1% Cu 0.5%, 
and are 600 A (26 fim) and 3000 A (10 /im) thick (wide) 
respectively. The junctions were patterned to a rectan- 
gular shape with area ^4 = 4x1 itm 2 by a self-aligned 
microfabrication process. 



We also define the resistance shift (8) in each cycle: 

-Rfinai -^initial 



8 = 



(^initial + -Rfinai) /2 



(C) 



IV. EXPERIMENTAL RESULTS 

The studied tunnel junction had an initial electri- 
cal resistance R = 57.8 fl and a resistance area product 
R x A = 230 fiyum 2 . No magnetoresistance was observed 
in our tunnel junctions, due to the loss of interfacial po- 
larization (20 A Ta layer deposited just below the bar- 
rier). In fact, the tunnel magnetoresistance of a TJ is 
known to exponentially decrease with the thickness of 
a non-magnetic layer inserted just below the insulating 
|- )arr j eB 28 i 29 anc j TMR then goes rapidly to zero within 
the first monolayers of the non-magnetic material. 

We measured CIS cycles with increasing J ma x, starting 
with a cycle up to imax = 30 mA (Fig. (2t; cycle start- 
ing at point S) giving CIS = 9.2% and 8 = -3.5%. No 
resistance switching was observed under the initial neg- 
ative current pulses (I p — — > — / ma x)- However, upon 
reversing the current one observes that for I p > 15 mA 
(where we define the positive critical current ; see Fig. 



EK) the resistance starts to decrease, a trend which be- 
comes increasingly enhanced (switching) with / p , up to 
I max = 30 mA. This switching is associated with elec- 
tromigration of metallic ions from the electrodes into the 
barrierjAiiii decreasing the effective barrier thickness and 
consequently the junction resistance. The previous ab- 
sence of R-switching under negative current pulses in- 
dicates an electromigration asymmetry with respect to 
the electrode/oxide interfaces, i.e. only ions from one 
such interface are actively participating in electromigra- 
tion. Physically such asymmetry arises not only from the 
different materials deposited just below (Ta) and above 
(CoFe) the insulating barrier, but also from the deposi- 
tion and oxidation processes during tunnel junction fab- 
rication. In particular the top electrode is deposited over 
an oxidized smooth surface, while a much more irreg- 
ular bottom electrode/oxide interface is experimentally 
observed^ Since the migration of ions into and out of 
the barrier should occur preferentially in nanoconstric- 
tions (higher electrical fields), one concludes that such 
ions likely belong to the Ta bottom electrode. The cur- 
rent density and electrical field at R-switching can be es- 
timated as j c ~ 0.375 x 10 6 A/cm 2 and E c ~ 3 MV/cm, 
respectively. 

Returning to Fig. |2k, the subsequent decrease of I p 
from +ima X to zero hardly affects the low resistance 
state. However, for I p < —15 mA (where we define the 
negative critical current /"), the resistance gradually in- 
creases until I p = — /max, recovering a significant frac- 
tion of the previous R-switching near +I max . This in- 
dicates that, under a reversed electrical field, many ions 
return to their initial sites. The subsequent change of 
I p from — /max to zero (to close the CIS cycle at point 
F) produces no significant change in resistance. How- 
ever, the final resistance mismatch (Rn na i < -Rinitiai; 
5 = —3.5%) indicates some irreversible effects in this CIS 
cycle (/max = 30 mA), associated with barrier degrada- 
tion. 

The voltage across the junction was also measured for 
each applied current pulse (/ p ), providing the V(/ p ) char- 
acteristic depicted in Fig. |2t> (hollow circles). If one 
uses Simmons' models to fit this curve with adequate 
thin TJ barrier parameters^ (barrier thickness t = 9 A, 
barrier height <j> = 1 eV), the quality of the fit is poor 
(dashed line in Fig. 0d), with large discrepancies near 
±/max- Also, the use of the Brinkman model for asym- 
metric tunnel junctions^! does not yield good fits. Such 
discrepancies near ± I max are related to localized heating 
inside the tunnel junction, as discussed below. 

We then performed CIS cycles with increasing 7 max , 
from 30 to 80 mA, in A/ max = 5 mA steps as shown 
for representative cycles in Fig. Notice the enhanced 
R-switching and R-recovering stages (versus / p ), occur- 
ring from /+ to / max and from I~ to — / m ax respectively. 
From these data one can obtain the CIS and (5-shift in 
each cycle, obtaining the corresponding dependence on 
/max as depicted in Fig. 0] The CIS coefficient rises with 
/max until ~65 mA (CIS = 57.4%), saturating for higher 
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FIG. 2: a) Current Induced Switching cycle for / max = 30 
mA, starting at point S and finishing at F. After each current 
pulse I p , the electrical resistance of the tunnel junction is 
measured under a low bias current, enabling us to obtain the 
depicted R(/ P ) cycle. Effective switching occurs between 
and +/ m ax, and resistance recovery between 1^ and — /max. 
b) Corresponding experimental (hollow circles) and simulated 
(dashed line) V(/ p ) characteristic. While applying the current 
pulse I p , the voltage across the tunnel junction is measured 
and a V(/ p ) characteristic obtained. 



FIG. 3: Selected CIS cycles performed with / max up to 80 mA. 
Notice the enhanced R-switching occurring under increasing 

/max • 



current pulses. On the other hand, S remains fairly small 
below /max ~ 60 mA (-0.4%), but increases rapidly for 
higher / max (S = -9.6% for Z max = 80 mA). The CIS 
increase with / ma x indicates that electromigrated ions 
are further pushed into the barrier (further lowering R) 
or/and more ions participate in the EM processes. Ul- 
timately irreversible damage occurs in the barrier, as 
reflected in the (5-shift enhancement for / max > 60 mA 
(Fig. 01 . 
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FIG. 4: Current Induced Switching coefficient and <5-shift as 
a function of maximum applied current. Large 5-shift values 
occur for / max > 60 mA, indicating progressive barrier degra- 
dation. 



FIG. 5: Resistance switching directions for a) 
AFM/FM/NM/I/FM (Mnlr/CoFe/Ta/AlO^/CoFe/NiFe) 
and b) AFM/FM/I/FM (Mnlr/CoFe/AlOx/CoFe/NiFe)^ 
tunnel junctions. 



V. DISCUSSION 

The observed resistance switching (R decrease) oc- 
curs only for positive current pulses in the here studied 
FM/NM/I/FM tunnel junctions (R-recovery occurs un- 
der negative I p ; see Figs. Et, EH and EJ , whereas in the 
previously studied FM/I/FM1& tunnel junctions switch- 
ing (recovery) occurs under negative (positive) currents 
(Fig. HJd; A = 2 x 1 fim 2 ). To explain such different be- 
havior one will compare electromigration direct and wind 
forces in Ta (NM) and CoFe (FM) layers. Using eq. 
we obtain: 

Z w (Ta) _ p(CoFe)v F (Ta) a t r(Ta) 
Z w (CoFe) ~ p(Ta)v F (CoFe) <7 tr {CoFe) { ' 

where vf is the Fermi velocity. Inserting the param- 
eters given in Table ^ 5 i 26 i 32 i 33 one obtains Z w (Ta) ~ 
0.07 Z w (CoFe). The wind force is then much larger in 
CoFe than in Ta layers and likely dominates electromigra- 
tion in the CoFe layers. On the contrary, because Ta is in 
an amorphous state (notice its high electrical resistivity 
in TableQ}, one expects the small electron mean free path 
to prevent large momentum gains by electrons between 
consecutive collisions. Using the value estimated previ- 
ously for Z w (Fe), one finds Z w (Ta) ~ —1.4 (rs Zd(Ta)). 
Remembering that the magnitude of the direct force is 
enhanced relatively to the wind force in nanoconstrictions 
(cq. see also below) and that the ballistic model over- 
estimates Z w , one expects the direct force to dominate 
in Ta. Thus, the likely cause for the observed difference 
in the R-switching directions is related to the dominance 
of different electromigration forces in Ta and CoFe. Con- 
firming this conclusion, tunnel junctions with Ta layers 
deposited just below and just above the insulating bar- 
rier (FM/NM/I/NM/FM; not shown) display the same 
current switching direction as those with only one Ta 



p (pQcm)^ a tr (A 2 )^ v F (cm/s 



CoFe 


17.1 


~3 


~2 


Ta 


154.0 


~6 


0.67 



TABLE I: Electrical resistivity, electron transport cross sec- 
tion for scattering and Fermi velocity used to estimate 
Z w (Ta)/Z w (CoFe). 



layer below the insulating barrier (FM/NM/I/FM). On 
the other hand, when a single NM Ta layer is deposited 
just above the barrier (FM/I/NM/FM), the R-switching 
direction is that of FM/I/FM tunnel junctions. 

Figure U3 (left scale) shows the CIS R(/ p )-cycle 
obtained at room temperature, with I max = 80 mA 
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FIG. 6: CIS cycle and corresponding V(ip) characteristic for 
1 max = 80 mA. Notice the decrease of |V| near zh 1 max ■ The 
dashed line depicts a V(I) curve calculated using Simmons' 
model. 



(CIS = 55.5%; 8 = -9.6%). Notice the R(/ p )-switching 
from 7+ = 15 mA to 7 max = 80 mA and resistance re- 
covery from I~ = —35 mA to -7 max = —80 mA. The 
V(7 p ) characteristic is also displayed (hollow circles; 
right scale), showing an anomalous plateau with a slight 
dV/d7 p negative slope for |7 p | > 30 mA. This effect can- 
not be explained by tunnel transport theories and is here 
related to heating inside the tunnel junction. Using our 
temperature dependent R-data^ the temperature inside 
the tunnel junction is estimated as ~600 K. Such high 
temperatures have also been observed in similar mea- 
surements performed in FM/I/FM tunnel junctions^ 

Heat generation in tunnel junctions arises from two 
processes^ usual Joule heating in the metallic layers and 
inelastic electron scattering upon ballistic tunneling. The 
steady-state heat equation can then be written ast^ti 
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where K is the heat conductivity, T is the temperature, x 
is the stack position, j = V/(RA) is the current density, 
V is the bias voltage and k n is the inelastic scattering 
electron mean free path. We obtained numerical results 
assuming that the current density is constant through- 
out the junction stack. The temperature at the bottom 
and top of the tunnel junction stack is assumed fixed at 
300 K. 

Our numerical results (Fig. |2j) indicate that large heat- 
ing can occur near the insulating barrier for high current 
densities. However, the temperature increase expected 
from the uniform case, j c = I c /A ~ 0.375 x 10 6 A/cm 2 
is negligible (~1 K; inset of Fig. EJ), and to reach 
600 K one needs j es t ~ 16 x 10 6 A/cm 2 . This corre- 
sponds to an effective area through which current flows 
A e ff — I c /jest ~ 0.1 /im 2 , i.e., about 2.5% of the total 
tunnel junction area. These results then suggest that 
j c is only an average value and that nanoconstrictions 
where the insulating barrier is thinner concentrate most 
of the current flowing through the junction. Such hot- 
spots have been observed in similar TJs by atomic force 
microscopy^ 

One can now understand the observed electromigration 
driven resistance changes in thin FM/NM/I/FM tunnel 
junctions with NM=Ta (amorphous; Fig. QJi). Under 
increasing positive current pulses (directed from the bot- 
tom to the top lead), the dominating EM direct force 
induced by the electrical field pushes Ta atoms into the 
barrier, a process thermally assisted by heating gener- 
ated by the high current densities flowing in nanocon- 
strictions. This rises the probability that an atom sur- 
mounts the energy barrier for migration Eb (see Fig. QJ, 
greatly enhancing atomic mobility. One notices that even 
a small barrier weakening (due to such migration) would 
considerably lower the tunnel resistance due to its ex- 
ponential dependence on barrier thickness 4& Using the 
Simmons' model we can calculate the resistance varia- 
tion due to a small barrier thickness reduction from t to 



FIG. 7: Simulation of heating processes inside the studied 
tunnel junction, under different electrical current densities 
(MA/cm 2 ). Inset: temperature increase as a function of cur- 
rent density passing through the junction. The lines show the 
current density needed for the temperature inside the junction 
to reach 600 K (1 MA= 10 6 A). 



t - St (St < t): 

R(t) - R(t - St) 



R(t) 

-^initial ~ ^half 
T^initial 



-B(<f>)8t 



(9) 



B((j>)St 



where B(<f>) = Q.72^/Jj2. For a CIS coefficient of ~ 60% 
one obtains a barrier thickness decrease St ~ 0.8 A. We 
can now plot the magnitude of the expected St decrease as 
a function of the maximum applied current 7 max (Fig. |H1 
using the experimental Timitiai and i?haif values), which 
follows the same trend as the CIS coefficient (Fig. EJ. 
In particular, a non-linear behavior (apparently expo- 
nential, as more clearly visible at low temperatures^) 
is observed for 7 max < 60 mA, that is, while the (5-shift 
is small and electromigration is mainly reversible. In 
atomic diffusion processes one often has2& Jp- oc F (x 
the position and t' the time). Therefore, in electromi- 
gration St oc ESt' , i.e. the barrier thickness decrease is 
proportional to the applied electrical field density and to 
the migration time St 1 . Following this simple analysis, 
one has (R(t) - R(t - St)) / R(t) oc E. The CIS effect 
then depends on how local electrical fields behave near 
nanoconstrictions and on its dependence on nanostruc- 
tural atomic rearrangements. 

Time dependent measurements (over 4 h) revealed 
that R remains practically constant both in its high and 
low state (not shown). This indicates that under a re- 
duced driving force, displaced Ta ions remain trapped 
in deep enough local energy minima inside lattice poten- 
tial barriers (E b 3> k B T), so that thermal fluctuations 
cannot return them to the electrodes. For example, in 



7 




- 0.00 

30 45 l max ( mA ) 60 75 90 



FIG. 8: Dependence of the effective barrier thickness decrease 
(St) on the maximum applied current pulse, as obtained from 
the CTS(I max ) curve (Fig. EJl and eq. El (for (j> = 1 eV). Inset: 
energy barrier for migration from Ta into the barrier (E^) 
and vice- versa (E^), in the first (left) and last (right) CIS 
cycles. 

the CIS cycle of Fig. QJi one observes that the low resis- 
tance state persists for I p current pulses from f=a + J max 
down to « However, when I p < I~ the driving 

force gets strong enough to return displaced ions back 
into their initial positions in the NM layer. However, the 
final resistance does not exactly reaches its initial value, 
indicating progressive barrier degradation. Such degra- 
dation should result from metallic ions that remain in the 
barrier after the CIS cycle is completed. We also notice 
that in the initial CIS cycle with I max = 30 mA (Fig. Et) 
one has 1+ « |-^~|- This indicates that the driving force 
for electromigration into and out of the insulating barrier 
is approximately equal, i.e. the lattice sites where ions 
migrate to are energetically similar. Furthermore, Fig. 
El (see dashed line) shows that J+ w 15 mA throughout 
all the CIS cycles performed, indicating that cycling does 
not alter the EM force inducing atomic migration from 
Ta into the barrier. In other words, the energy barrier 
which the Ta ions surmount when migrating into the bar- 
rier is kept constant (inset of Fig. |HJ) . This contrasts with 
electromigration in the opposite direction, where \I~\ in- 
creases with cycling (Fig. see dotted line). The force 
needed to return ions back has to be increased (inset of 
Fig. |HJ , indicating that Ta ions migrating under increas- 



ingly higher current pulses are pushed further inside the 
barrier, and are thus more difficult to return to the elec- 
trode. 



VI. CONCLUSIONS 

We studied the Current Induced Switching effect on 
low resistance (7 A barrier) CoFe/Ta/A10a;/CoFe tun- 
nel junctions. The CIS coefficient increased with increas- 
ing maximum applied current pulses, reaching ~60% for 
I max = 80 mA. Such effect is controlled by nanostruc- 
tural rearrangements at the electrodes/barrier interfaces, 
due to ion electromigration (reversible and irreversible). 
When high currents are applied, one observes large irre- 
versible resistance decreases. The V(/ p ) characteristics 
showed an anomalous behavior when | I max | > 65 mA 
due to heating effects inside the tunnel junction, show- 
ing that the CIS effect is thermally assisted. The analysis 
of these effects shows that nanoconstrictions indeed con- 
centrate most of the tunneling current through the bar- 
rier, forming local hot-spots. One further demonstrates 
that the R-switching direction is related to a competi- 
tion between the electromigration contributions due to 
direct and wind forces: the direct force dominates elec- 
tromigration in Ta layers, whereas the wind contribution 
is dominant in CoFe. 

Finally, please notice that, although the results pre- 
sented here concern a single FM/NM/I/FM tunnel junc- 
tion, they are reproduced when measuring other TJs from 
the same deposition batch. Particularly, the dependence 
of the CIS coefficient on maximum applied electrical cur- 
rent is quite similar in different tunnel junctions. The 
current switching direction is always the same for the 
same TJ-structure. 
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